Complete integrability of the coupled KdV— mKdV system 
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Abstract. The coupled KdV— mKdV system arises as the classical part of 
one of superextensions of the KdV equation. For this system, we prove its 
complete integrability, i.e., existence of a recursion operator and of infinite 
series of symmetries. 



Introduction 

There are several supersymmetric extensions of the classical Korteweg-de Vries 
equation (KdV) || g). One of them is of the form (the so-called N = 2, A = 1 
extension [Q]) 

U t = —U3 + 6u7il — 3(fi(fi2 — 3lp1p2 ~ 3llW3 — ZW\W2 + 3uiw 2 + 6uwwi 

+ 6ip(p\w — 6tpipiw — 6(pijjwi, 

(ft = —^3 + 3ipui + 3tpiu — 3-02 w — 3-0iwi + Zipiw 2 + &(fww\, 

ipt — —03 + 30ui + 30im + 3(f2W + 3(fiwi + 3ij)iw 2 + 6ipwwi, 

w t = — W3 + 3w 2 wi + 3uwi + 3uiw, 

where u and w are classical (even) independent variables while <p and ip are odd 
ones (here and below the numerical subscript at an unknown variable denotes it 
derivative over x of the corresponding order). Being completely integrable itself, 
this system gives rise to an interesting system of even equations 

Ut = — W3 + 6uui — 3uiws — 3wiW2 + 3uiw 2 + 6uwwi, 

w t = —w 3 + 3w 2 wi + 3uw\ + 3ui%u, (1) 

which can be considered as a sort coupling between the KdV (with respect to u) 
and the modified KdV (with respect to w) equations. In fact, setting w = 0, we 
obtain 

Ut = — U3 + 6uu%, 

while for u = we have 

Wt — — W3 + 3w 2 wi. 
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In what follows, we prove complete integrability, cf. j2j, of system ([!]) by estab- 
lishing existence of infinite series of symmetries and/or conservation laws. Toward 
this end we construct a recursion operator using the techniques of deformation 
theory introduced in Q and extensively described and exemplified in |Q . 

In the first section of the paper the theoretical background is introduced. The 
second section deals with particular computations and description of basic results. 

1. Geometrical and algebraic background 

Here we briefly describe the geometrical theory of partial differential equations 
and algebraic foundations of computational approach to recursion operators 
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Let 7r: E — > M be a locally trivial vector bundle and ir^ : J k (^) — > M, k = 
0, f , . . . , oo, be the bundles of its fc-jets. A (nonlinear) partial differential equation 
(PDE) of order k is a submanifold £ C J k {^), k < oo. Its Ith prolongation is a 
subset £ l C J k+l (n). There exist natural mappings nk+l+i,k+l '■ ~~ > £\ an d £ 
is said to be formally integrable, if all £ l are smooth manifolds while iik+i+i,k+i are 
smooth fiber bundles. Below, only formally integrable equations are considered. 

The inverse limit £°° C J°°(7r) of the system {£ , irk+i,k+l-i} is called the 
infinite prolongation of £ and we consider the bundle n£ : £°° —> M . This bundle 
enjoys the following characteristic property: 

Proposition I (see ||). Let TVi\ E t — » M, i — 1,2, be two locally trivial vec- 
tor bundles and A: r(?ri) — » r(7T2) be a linear differential operator acting from 
sections of m to sections of 7T2 ■ Then there exists a unique differential operator 
CA: r(7r£(7ri)) -> r(7r|(7r 2 )) such that 

j 0O (s)*oCA = Aoj 0O ( s )* (2) 

for any formal solution s of the equation £. The correspondence A i— * CA is 
C 00 (f °) -linear and complies with the composition of differential operators: 

C(A a oAi) = CA a oCAi, (3) 

where A±: T(tti) — > r(?r2), A2 : r(7T2) — > r(7T3) are linear differential operators, 
773 : E3 — > M being a third vector bundle. 

As a corollary of Proposition |l|, we get 

Proposition 2. The bundle its possesses a natural flat connection. 

PROOF. It suffices to take the trivial bundle 1m- M x R —> M for the bun- 
dles 7Ti and 7T2 and an arbitrary vector field X for the operator A. Flatness is a 
consequence of (||). □ 

Let us denote by D(iV) the C°°(iV)-:module of vector fields on a manifold N. 

Definition 1. The connection C: D(M) — ► D(£°°) is called the Cartan con- 
nection on £ 00 . 

Denote by CD(£°°) C D(£°°) the horizontal distribution on £°° with respect to 
the Cartan connection (the Cartan distribution) and by De(£°°) the normalizer of 
CD(£°°) in D(£°°). Then, since C is flat, CD{£°°) is integrable in a formal Frobenius 
sense and thus CD(£°°) is an ideal in Dc(£°°). 

Definition 2. The quotient Lie algebra sym£ = D C (£ 00 )/CD(£ 00 ) is called 
the algebra of higher symmetries of the equation £. 
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The Cartan connection in tt £ determines the splitting 

D(£°°) = D v (£°°) © CD(£°°), (4) 

where D v (£ co ) denotes the module of 7r£-vertical vector fields, and for any coset S £ 
sym £ there exists a unique vertical representative. In its turn, such a representative 
is uniquely determined by a section ip £ T(ir£(ir)) (generating section) satisfying 
the defining equation 

ts<P = (5) 

and vice versa. Here £ £ is the universal linearization operator for £ restricted to 
£°° . Due to this fact, we shall identify the solutions of (|) with higher symmetries 
of£. 

The connection form U £ (the structural element of £) of the Cartan connection 
C is an element of the module D"(A 1 (f 00 )) of A 1 (£' 00 )-valued vertical derivations. 
Thus, we can introduce an operator d £ : D v (A i (£°°)) -> D"(A i+1 (£°°)) defined by 

d £ Q=lU £ ,Ql OeD"(A ! (r)), 

where [■, ■] is the Frolicher-Nijenhuis bracket. Since the Cartan connection is flat, 
one has \U £ , U £ \ — 0, from where it follows that d £ o d £ = 0. Thus we obtain 
a complex (D v (A l (£°°)),d £ ), whose cohomology is called the C-cohomology of £ 
and is denoted by H'(£). It is easy to see that Hq(£) — sym£, while H^{£), by 
standard reasons, is identified with classes of nontrivial infinitesimal deformations 
of U £ (or, which is the same, of the equation structure). 

If VL and 6 are elements of B v (A l {£°°)) and D V (A J (£°°)) respectively, their 
contraction £1 _i 9 is defined as an element of (D V (A' ,+: >~ 1 (£ 00 )) This operation is 
inherited by the C-cohomology groups. In particular, if p £ sym£ and 1Z £ Hq(£), 
then tp _ilZ — IZtp is a symmetry again. In other words, the module H^{£) acts on 
the Lie algebra of higher symmetries. 

Let CA 1 (£°°) C A 1 (f°°) be the submodule of one-forms on £°° vanishing on 
the Cartan distribution. Then one has the direct sum decomposition 

A x (f °°) = CA 1 ^ 00 ) e A 1 ^ 00 ) 

dual to (^), where A^(£°°) is the submodule of horizontal forms. This splitting is 
also inherited by H^(£°°) and an element 1Z £ H^(£°°) acts nontrivially on sym£ if 
only it corresponds to a derivation from D U (CA 1 (£°°)). Moreover, it can be shown 
that \md £ R D'"(CA 1 (£ 00 )) = and consequently nontrivial actions can be found 
by solving the equation 

d £ (R) = 0. (6) 

Solutions of (|^) are called recursion operators for symmetries. Any recursion oper- 
ator TZ is uniquely determined by an element ton £ CA 1 (^°°) <E> T(ir £ (ir)) satisfying 
the defining equation 

4 X W)=0, (7) 

where is the extension of the operator i £ to the module CA X (£ °°) ®V{^k £ {^k)). If 
p £ r(Tr £ (Tr)) is a symmetry and S v £ D v (£°°) is the corresponding vertical vector 
field, then the action of TZ on p is given by 

TZp = S v -iLu n . (8) 
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Local coordinates. Let U C M be a coordinate neighborhood in M such that 
the bundle 7r trivializes over 14, x%, . . . ,x n be local coordinates in Li and u , . . . , it™ 
be coordinates along the fiber in a given trivialization. Then the adapted coordinates 
in ir^-(U) C J°°{n) are given by the functions u 3 a , j = l,...,m, <r = ii . . - ijfe, 
1 5- i« < ^, uniquely defined by 



joo(/r(«i) = 



. . . ctoi 



for any local section / = (J 1 , . . . , J" 1 ) e r(; 
J°°{tt) is given by 



Then the Cartan connection in 



(9) 



where Di are the so-called total derivatives. The structural element in this case is 
given by the formula 

d 



dui 



(10) 



where LJ J a = du J a — c?Xj are the Cartan forms constituting a basis in the 

module CA 1 ( J°°(ir)). The forms u>l may be rewritten as u>%. — dcu 3 a , where dc is 
the Cartan differential acting on / S C 00 (J 00 (7r)) by 

This differential restricts to any submanifold £°° C J°°(tt) and for infinite prolon- 
gations ( |To| ) transforms to 



_9_ 

9u/ 



(11) 



where {it/} spans the set of internal coordinates in f 00 . 

If an equation £ C J k (tt) is given by the system of equalities 

F (x%) . . . , x n , . . . , w^., . . . ) 0, 

< 

F r (xi, x n ,..., <,...)= 0, 

then its infinite prolongation is defined by 

D a F a = 0, a = l,...,r, 0< |ct| < oo, 

where = o- • -oDi k for a = i\ . . .i^. Then the universal linearization operator 
corresponding to this system is of the form 



(12) 



The operator is is obtained from ([L2h by rewriting it in internal coordinates. 
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Example (Evolutionary 1 + 1 equations). Consider a system £ of evolution- 
ary equations 



'v} = F 1 (x,t,v. 1 ,...,vF), 



(13) 



uf = F r (x,t,u 1 



,...,u k ), 



where x = Xi, t = %% and u\ — d v? /dx\ Then the functions x, t, . . . , u\ , ■ ■ ■ can 
be taken for internal coordinates on £°° . The total derivatives are written down as 



1,3 l 1,3 

in these coordinates. The Cartan forms on £°° are 

>jj\ = du\ - u\ +1 dx - D l x {Fi) dt, 
while the structural element for £ is given by 

d 



du J , 



(14) 



1,3 1 

The restriction 1% of the universal linearization operator to £°° is of the form 



Is 



y®El D i x -vD t 



where E is the identity matrix. So, a vector function ip = (cp ,. . . ,cp m ), (p 3 = 
Lp 3 (x, t, . . . , uj , . . . ) , is a symmetry of £ if and only if 

f)TP a 

Y^—D 1 ^ = D t p a (15) 

1,3 dU l 

for all a = 1, . . . , m. The corresponding vector field is given by the formula 

~3 9U 1 

In a similar way, recursion operators are determined by vector-valued forms 
uj n = (w^, . . . , uj%), where uj 3 r = J2t a ipa^f, i/4t e C 00 ^ 00 ), satisfying the equa- 
tions 



1,3 Ul 

To compute the left and right sides of ([l6]), it suffices to note that 
D m <4 = wf +1) D t u4 = D l x d c Fi =D l x J2 q^»%- 

OL.be $ 



(16) 



Remark 1. Let lu = (w 1 , . . . , uj" 1 ), uj j = J2i a ''Pa 10 ?' be a vector-valued Car- 
tan form on £°° . Then for any vector- valued function ip — (p 1 , . . . , ip m ) on £°° the 
action 72- w : <p <—>■ TZujip is defined by IZ^p — 9 v -ilo. In local coordinates, this action 
is expressed by the formula 
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Operators of this type (i.e., expressed in terms of total derivatives) are called C- 
differential (or total differential) operators. 

Remark 2. As it was mentioned above, operators of the form ([l7]), provided 
lo satisfies (|l6|), take symmetries of the equation at hand to symmetries of the same 
equation. In other words, one has TZ^ : ker£g — > ker^g. Under not very restrictive 
conditions on the equation £, this is equivalent to the operator equality 

ieoH u = Aoi e , 

where A is a C-diffcrential operator. Taking formally adjoint, one obtains 

Ktat £ =t £ aA*, 

which means, that A* is a recursion operator for generating functions of conserva- 
tion laws [fit 



Nonlocal setting. In practice, when solving (Q), one usually finds no nontriv- 
ial solutions, though the equation £ may possess recursion operators. 

Example (The Burgers equation) . Consider the equation 

U t = U xx + Wi- 
lt is known to possess a recursion operator of the form 

K = D x + -u + -u 1 D- 1 . (18) 

Nevertheless, the only solution of the equation 

t [ £ x] LO = {Dl + u D x + ui - D t )w = (19) 

for lo = Vwo + • ■ • ipk^k is cvloq, agl, which provides the trivial action 1Z U : ip \— > 
aip. 

To resolve this apparent contradiction, let us extend the algebra C 00 ^ 00 ) with 
an additional element U-\ and set 

D x U-i = M , 

1 2 

D t u-i = ui + -u , 

dcU—i = LO—i = — uq dx — (ui + — dt. (20) 

Then, solving ( |l9| ) for to = + ipo^o + • • ■ V'fc^fci we obtain a two-parametric 

solution 

LO = OtLOQ + (3ft, ft = LOl + —UqLOo + — 

Then the action ip i— > TZ^ip — 9 V _i Q coincides exactly with (|T^). 

This example reflects a general scheme of computations which arises in a lot 
of applications and is used in Section |[ Namely, in search of recursion op- 
erators for symmetries we extend the algebra C°°(£°°) with a new set of vari- 
ables w , . . . , w r , . . . (so-called nonlocal variables) and respectively extend the total 
derivatives to vector fields 

A = A + $>^> i = l,...,n, (21) 
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in such a way that 

[A , Dj] = [A , Xj] + [X t , Dj] + [X t ,X,]= 0, (22) 

where Xi = X?d/dw a . 

Having a solution Xx, . . . ,X n of (pl|), we obtain an integrable distribution on 
the space £°° x R^, where N is the number of nonlocal variables (the case N = oo 
is included). The projection r: £°° = £°° x M. N — > £°° is called a covering over 
£ and iV is called its dimension (for an invariant geometrical definition s ee Q ). 
Similar to the local case, we define the Lie algebra sym T £ = Dc(£ 00 )/CD(£°°) of 
nonlocal r-symmetries. 

We introduce Cartan forms 

a 

0i = dw j -^XfdXi, j 1 V. 

corresponding to nonlocal variables on The module of all Cartan forms on 

£°° is denoted by CA 1 (£°°). We also extend the universal linearization operator £g 
to £°° just by changing the total derivatives A to A- Let us now consider two 
equations, associated to this extension: 

e £ tp = o and = o, 

where (p e T((t: £ o t)*tt) and Q e r((7r f o t)*tt) ® CA 1 (f°°). Solutions of the first 
equation are called r-shadows of nonlocal symmetries, while solutions of the second 
one are said to be r-shadows of recursion operators in the covering t. The following 
result establishes relations of shadows to symmetries and recursion operators: 

Theorem 1 (sec @, §). Let t:£~^ -> 6e a covering. Then: 

1. If ip is a T-shadow, then there exists a covering 

f:W->£~°°^£°° 

such that Lp reconstructs up to a nonlocal f -symmetry. 

2. If Lp is a nonlocal r-symmetry and 1Z is a recursion operator shadow, then 
IZip is a symmetry shadow. 



Remark 3. Among all one-dimensional coverings over a given equation there 
exists a special class consisting of those ones, for which the fields X\, . . . , X n in (22) 
are independent of nonlocal variables (so-called abelian coverings) . To any such a 
covering, one can put into correspondence a differential form on £°°: 



1=1 



This form is closed with respect to the so-called horizontal differential dh = Cd 
(cf . Proposition 0) . Vice versa, to any such a form there corresponds a covering of 
the above mentioned type. Moreover, two closed forms determine the same class 
in the cohomology group H 1 (dh) if and only if the corresponding coverings are 
equivalent. In particular, if n = 2, the group H 1 (df l ) coincides with the group 
of conservation laws of the equation £ |fTof| . Thus, to construct a covering under 
consideration is the same as to find a conservation law. This fact is used in the 
computations below. 
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2. Basic computations and results 

In this section we shall discuss the complete integrability of the KdV-mKdV 
system given in (|l]), i.e., 

u t — — U3 + &uu\ — iww?, — iw\W2 + 3u\w 2 + Guwwi, 

w t — — 1^3 + 3w 2 w\ + 3uwi + 3uiw. (23) 

In order to demonstrate the complete integrability of this system, we shall construct 
the recursion operator for symmetries of this coupled system, leading to infinite 
hierarchies of symmetries and, most probably, of conservation laws. Due to the 
very special form of the final results, it seems that integrability of this system, 
which looks quite ordinary, has not been discussed before or elsewhere. In order to 
do this, we shall discuss conservation laws in Subsection 2.1 leading to the necessary 
nonlocal variables. 

In Subsection |2.2| w e shall discuss local and nonlocal symmetries of the system, 
while in Subsection |2.3| we construct the recursion operator or deformation of the 
equation structure (|14|). 

2.1. Conservation laws and nonlocal variables. Here we shall construct 
conservation laws for ( ^3| ) in order to arrive at an abelian covering of the coupled 
KdV-mKdV system as was shown for Burgers equation (poj). So we construct 
X = X(x, t,u, . . . , w . . . ), T = T(x, t,u, . . . , w . . . ) such that 

D x {T) = D t (X) (24) 

and in a similar way we construct nonlocal conservation laws by the requirement 

D x {f)=D t (X), (25) 

where D* is defined by (||); moreover X, T are dependent on local variables x, 
t, u, . . . , w, . . . as well as the already determined nonlocal variables, denoted here 
by or p t! », which are associated to the conservation laws (X, T) by the formal 
definition 

D x {p*) = (p*) x = X, 
D t {p*) = (p*)t = T. 
Proceeding in this way, we obtained the following set of nonlocal variables 

P0,1> P0,2, Pi, Pl,D Pl,2> P2,l) P3, P3,l, P3,2, P4,l, P5, (26) 

where their defining equations are given by 
(Pl)x = u, 

(Pi)t = 3it 2 + 3uw 2 — ii2 — 3wiV2, 

(P0,l)x = W, 

(po,i)t = 3uw + w 3 -w 2l 

(P0,2)x =Pl, 

(P0,2)t = ~6p 3 - Ui, 

{Pi.i)x = cos(2p ,i)Pi w + sin(2p ,i) w2 i 

(pi,i)t = cos (2po,i)(3piwif + Piw 3 — p\W2 + uwi — u\w — w 2 W\) 
+ sin(2p ,i)(4MW 2 + w 4 — 2ww 2 + w\), 
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{Pi,2)x = cos(2p ,i)w 2 - sin(2p . 1 )p 1 w, 

(Pi,2)t = cos(2p ,i)(4uu; 2 + w 4 - 2ww 2 + w\) 

+ sin(2po,i)(— Spiuw — piw 3 + piw 2 — uwi + Uiw + w 2 wi), 

(P2,i)x = (4cos(2p 04 )p M w 2 -4sin(2p 0!l )p 1 p 14 u; + w7(p 2 - 2u + w 2 ))/2, 

(P2,i)i = (4cos(2p ,i)Pia(4Mw 2 + w 4 - 2ww 2 + wf) 

+ 4sin(2po,i)pi,i(— 3piuw — piw 3 + P1W2 — »i + uiiu + ui 2 wi) 
+ 3p 2 uw + p\w 3 — p\w2 + 2piuw\ — 2pimw — 2piw 2 wi — 8u 2 w 

- uw 3 + 2uw 2 — 2u\Wi + 2u 2 w + w 5 + 2>w 2 w 2 )/2, 
(?3)x = (-u 2 - uw 2 + ww 2 )/2, 

{ps)t — (—4m 3 — 9w 2 w 2 + 2uu 2 — 2>uw 4 + lluww 2 — uw 2 — u\ + u\ww\ 

+ 4u 2 w + 6w W2 + 3w w 1 — WW4. + — w 2 )/2, 

(P3,i)x = (cos(2p ,i) l0 (Pi ~ 6 Pi u + 39piw 2 - 24pi ; ipi i2 w + 12p 3 + 6ui) 

+ 2 sin(2p ,i )w(12pipi,ipi,2 + 18piWi + 2uj 3 + 3u> 2 ) 

+ 6pi, 2 w(-p 2 + 2u - w 2 ))/12, 
(P3,2)x = (2cos(2p ,i)w(12piPi,iPi,2 - 18piu>i - 2w 3 - 3w 2 ) 

+ sin(2p ,i)w>(p 3 - 6piu + 39piw 2 + 24pi i ipi ;2 u> + 12p 3 + 6m) 

+ 6 Phl w(-pj + 2u - w 2 ))/l2, 
(P4,i) x = (8cos(2p ,i)w(pi"Pi,2 + 12pip 2 1 pi i2 - Qpipi, 2 u + 3p!p lj2 w 2 

- I2pi t ip1 2 w + l&pi^uw - 4pi j iw 3 - 6pi,iw 2 + 12pi,2P3 + 6pi, 2 ui) 
+ 8 sin(2p ,i) u; (PiPi,i + 1 2 PiPi,iPi,2 ~ 6pip M u + 3pip M u; 2 

+ 12p' 2 1 pi j2 w + 12pi i ip 3 + 6^1,1^1 - 18pi j2 uw + 4pi, 2 w 3 + 6^1,2^2) 
+ w(-pj - 24p 2 p 2 ;1 - 24p 2 lP l 2 + \2p\u - Qp\w 2 - 48 PlP3 

- 24pi«i + 48p 2 s u - 2Ap\ 1 w 2 + 48p 2 2 M - 2Ap\ a w 2 

- 60u 2 + Uuw 2 + 24u 2 - 13w 4 + 6uw 2 ))/48, 

iPb)x = (12m 3 + 24m 2 w 2 — 6uu 2 + 6uw 4 — 30uww 2 — 3ii2W 2 — 8u> 3 u>2 + 6wu)4,)/6. 

In the previous equations, we skipped explicit formulas for (p3s)t, (ps, 2 )t, 
(P4,i)t, and (p5)t, because they are too massive, though quite important for the 
setting to be well defined and in order to avoid ambiguities. The reader is referred 
to the Appendix for them. 

It is quite a striking result that functions cos(2p ,i), sin(2po,i) appear in the 
presentation of the conservation laws and their associated nonlocal variables. 

We should note that pi, po,i) P3, P5 arise from local conservation laws and 
we shall call pi, po,i, P3, P5 nonlocalities of first order. 

In a similar way we see that po,2, Pi, 2 arise from nonlocal conservation 
laws, where their x- and t-derivatives are dependent on the first order nonlocalities. 
For this reason p ,2, Pi,i, Pi, 2 are called nonlocalities of second order. 

Proceeding in this way p 2 ,i, P3,i, P3,2, P4,i constitute nonlocalities of third order. 
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2.2. Local and nonlocal symmetries. In this section we shall present re- 
sults for the construction of local and nonlocal symmetries of system (|3|). In order 
to construct these symmetries, we consider the system of partial differential equa- 
tions obtained by the infinite prolongation of ( p3| ) together with the covering by 
the nonlocal variables 

Pt),l, P0,2, Pi, Pl,l, Pi, 2, P2,l, P3, P3,l, P3,2, P4,l, P5- 

So, in the augmented setting governed by (p3[), their total derivatives and the 
equations given in Subsection 2.1 we construct symmetries Y — (Y U ,Y W ) which 
have to satisfy the symmetry condition 

i £ Y = 0. 

From this condition we obtained the following symmetries 

5o,l ; 51,1) 51,2) 51,3, 52,1, 53,1, Y 3 ,2, 53,3, 

where generating functions F"^ , Y™^ are given as 
Yq A = 3t{ t ouui + 6uwwi + 3uiiu 2 — — 3ww^ — 3wiW2) + xui + 2u, 
Yq°i — 3t(3uwi + 3uiw + 3w 2 wi — W3) + xw\ + w, 

Yi™ 1 =w 1 , 

Y™ 2 = cos(2p ,i)(2uw — w 2 ) + sin(2p ,i)(ui + 2ww\), 
Y™ 2 = - cos(2p ,i)w - sin(2p ,i)wi, 
Yi. 3 = cos(2po,i)(wi + 2wwi) + sin(2p ,i)(— 2mu + io 2 ), 
Y™ 3 = - cos(2p ,i)wi + sin(2p ,i)w, 

Y 2 U A = (2cos(2p ,i)(pi,iwi + 2p lt iwwx - 2pi <2 uw +^1,2^2) 
+ 2sin(2p ,i)(— 2pi,iuw + pi,iw 2 - pi, 2 u i - ^Pi,2 ww i) 
+ 2piuw — P1W2 + 2uwi + 3u\w + 2w 2 w\ — W3)/2, 

Y 2 W i = (2cos(2p ,i)(-pi4 w i +Pi,2u) +2sin(2p ,i)(pi,i« + Pi,2i0i) 

- piu + ui + wwi)/2, 

Y31 = (6uui + 6uwwi + 3uiw 2 — W3 — 3wu>3 — 3iyiio 2 )/3, 
5^"i = {3uwi + 3u\w + 3w 2 wi — ws)/3, 

Y 3 U 2 = {cos(2po.i)(—2pi uw + p\w 2 — ipiuwi — 6piUiw — Apiw 2 w\ + 2p\w-$ 
+ Bpi^pi^ui + 16pispi, 2 wwi - 8p\ 2 uw + Ap\ 2 w 2 — Ap 2 .iU\ — 8p 2y iwwi 
+ 10u 2 w + 6uw 3 — 8uw 2 — 14uiWi — 8u 2 w — llw 2 w 2 — lAww 2 + 2W4) 
+ 2sin(2p ,i)(-8pi,ipi,2uw + ^Pi,iPi,2W 2 - 2p\ 2 u 1 - Ap\ 2 wwi + 4p 2 ,iuw 

- 2p 2 ^\w 2 + 6uui + lOuwwi + 3u\w 2 — U3 + 2w 3 wi — 3ww3 — 5wiW 2 ) 
+ 4p>x )2 {2piuw — piw 2 + 2uw\ + 3uiw + 2w 2 wi — u>3))/8, 

Y 3 W 2 = (cOs(2poa)(Pl M - 2piUl - 2PXWW1 - 8^1,1^1,2^1 + 4p?,2 u + 4^2,1^1 

- 4m 2 - 3uw 2 + 2u 2 + 4ww 2 + 2w\) 

+ 2 sin(2j>o,i)(4pi,iPi.2W + 2p\ 2 w\ — 2p 2 \u — 3uw\ — 3uiw — 3w 2 w\ + W3) 
+ 4pi,2(-pi« + ui + wwi))/8, 
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y 3 " 3 = (2cos(2p ,i)(2p 1)1 ui + Ap x l ww\ — 4p 2 ,iuw + 2p 2yl w 2 - 6uui 

— Wuwwi — 3u±w 2 + 113 — 2w 3 wi + 3ww3 + 5wiw 2 ) 

+ sin(2po : i)(~2p 2 uw + p\w 2 ~ Apiuwi — &p\Uiw — Ap\W 2 w\ + 2piW3 

— 8p\ jura + Ap 2 x w 2 — 4p 2 ,iUi — 8p 2 \WW\ + 10u 2 w + 6uw 3 

— 8uw 2 — IAuiWi — 8u 2 w — llw 2 w 2 — lAww 2 + 2-1/74) 

+ 4 : pi i i(2piuw — pitu 2 + 2uw\ + 3uiw + 2w 2 wi — W3))/8, 
Y™ 3 = (2cos(2poj)(—2p 2 x wi + 2p 2: iu + 3uwi + 3uiw + 3w 2 wi — W3) 
+ sin(2p ,i)(Pi u _ 2piui — 2piww\ + Ap\ s u + 4^2,1^1 — 4u 2 

— 3uw 2 + 2u 2 + Aww 2 + 2w\) 
+ 4pi,i(-piu + ui + wwi))/8. 

2.3. Recursion operator. Here we present the recursion operator 1Z for sym- 
metries for this case obtained as a higher symmetry in the Cartan covering of system 
of equations (|l|) augmented by equations governing the nonlocal variables (^6| ) . As 
explained in the previous section, the recursion operator is in effect the a deforma- 
tion of the equation structure ( |ll|) . 

As demonstrated there, this deformation is a form-valued vector field (or a 
vectorfield-valued one-form) and has to satisfy 

lf ] TZ = 0. (27) 

In order to arrive at a nontrivial result as was explained for Burgers' equation too 
(c.f. Example [l]), we have to introduce nonlocal variables 

PQ,l, P0,2> Pi: Pl,l) V\, 2, P2.1, P3, P3.1, P3,2, P4,l, P5 

and their associated Cartan contact forms 

W P0,1' W PO,2J ^Pll W Pl,l' W Pl,2J UJ P2,1> W P3' W P3,1! W P3,2J W P4,1' W P5 ' 

The final result, which is dependent on the nonlocal Cartan forms 

k'po.l' W Pl ' ^Pl.l' W Pl,2> 

is given by 

d d 

ll = K u — + K w — + ..., (28) 
ou aw 

where the components R u , R w are given by 

Ru = (J U2 (— 1) + ^«(4w + w 2 ) + cj W2 (—2w) + uj Wi {~wi) + uj w (3uw - 2w 2 ) 
+ Wpi >2 (— cos(2p ,i)(ui + 2wwx) + sin(2p ,i)(2uw - w 2 j) 
+ uo Pl l (cos(2po : i)(-2mui + w 2 ) - sin(2p ,i)( u i + 2ww{)) 
+ ui Pl (2ui + wwi) + u! po 1 (2p\uw — p\w 2 + 2uwx + 3uiw + 2w 2 w\ ~ W3), 
R w = Wtu 2 (— 1) + u w (2u + w 2 ) + lo u (2w) 
+ oj Pi 2 (cos(2p ,i) w i - sin(2p ,i) u ) 
+ uj Pia (cos(2p ,i) u + sin(2p ,i)wi) 

-I- uj Pi (wi) + uj P01 {-p\U + ui+ wwi). (29) 
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We shall now present this result in a more conventional form which appeals 
to expressions using operators of the form D x and D~ x . In order to do this, we 
first split (^9|) into the so-called local part and nonlocal parts, consisting of terms 
associated to u) U2 , u> u , lj W2 , u) Wl , u> w and those associated to ui Pl 2 , ui Pl 1 , u) Pl , u> Pol 
respectively. The first part will account for D x presentation, while the second one 
accounts for the D~ x part. 

Due to the action of contraction 9 ip _i7i, the local part is given by the following 
matrix operator: 



-Dl + 4u 
2w 



-2wD 2 — w\D x + 3uw — 2w2 
-Dl + 2u + w 2 



The nonlocal part will be split into parts associated to u> Pl , uj Po 1 and oj p1 2 , ui Pl 1 , 
respectively. The first one is given as 

(2u± + ww^D^ 1 (2piuw — P1W2 + 2uw\ + 'Suiw + 2w 2 w\ — W3)D~ 1 
wiD^ 1 (—piu + ui+ww^D^ 1 

To deal with the last part, let us introduce the notation: 

A\ = cos(2po,i)(— 2mu> + W2) — sin(2po.i)(ui + 2wwi), 

A 2 = cos(2p ,i)w + sin(2p ,i)^i, 

Bi = — cos(2po,i)(wi + 2wwi) + sin(2po,i)(2uu) — W2), 

B 2 = cos(2p ,i)wi - sin(2p ,i)u, 

being the coefficients at u> Pl l and oj P12 in (p9|). 

According to the presentations of {jp\,i) x and (pi^x, i-e., 

{Pi,i)x = cos(2poa)Pi w + sin(2p ,i) w2 
(Pi,2)x = cos(2p ,i) w2 - sin(2p ,i)Pi w 7 
we introduce their partial derivatives with respect to po t i, pi, and w as 

ai = — 2p\w sin(2po,i) + 2w 2 cos(2j>o,i), 

a 2 = wcos(2p ,i), 

a 3 =pi cos(2p ,i) + 2u;sin(2p ,i), 

0\ = — 2w 2 sin(2po,i) — 2piW cos(2po,i), 

02 = -twsm(2p ,i), 

f3 3 = 2wcos(2p ,i) - pi sin(2p ,i)- 

From this we arrive in a straightforward way at the last nonlocal part of the recur- 
sion operator, i.e., 

AiD x 1 a 2 D x 1 AiD^iaxD- 1 +a 3 ) 
AzD-^D- 1 A 2 D x 1 {a 1 D x 1 + a 3 ) 

So, in the final form we obtain the recursion operator as 

~-Dl+4:U + w 2 -2wD 2 x - w x D x + 3uw - 2w 
2w -D 2 x + 2u + w 2 

{2u\ + wwijD^ 1 (2piuw — P1W2 + 2uw\ + 3uiw + 2w 2 w\ 
wiD x x (—piu + Ux + ww\)D~ 1 

A 1 D- 1 a 2 D- 1 AxD^iotxD- 1 + a 3 
A 2 D~ 1 a2D x 1 A 2 D x 1 (a 1 D x 1 + a 3 



n 



+ 



+ 



B 1 D x 1 p 2 D x 1 BxD-^PxD- 1 + fa) 
B 2 D x 1 faD x 1 B 2 D x \p 1 D x 1 + fa) 



COMPLETE INTEGRABILITY OF THE COUPLED KDV-MKDV SYSTEM 13 

+ [B 2 D-^ 2 D-^ B 2 D- 1 (0 1 D- 1 +(3 3 )_ ■ 

3. Conclusion 

We gave an outline of the theory of deformations of the equation structure of 
differential equations, leading to the construction of recursion operators for sym- 
metries of such equations. The extension of this theory to the nonlocal setting of 
differential equations is essential for getting nontrivial results. The theory has been 
applied to the construction of the recursion operator for symmetries for a coupled 
KdV-mKdV system, leading to a highly nonlocal result for this system. Moreover 
the appearance of nonpolynomial nonlocal terms in all results, e.g., conservation 
laws, symmetries and recursion operator is striking and reveals some unknown and 
intriguing underlying structure of the equations. 

Appendix 

Here we present explicit formulas for (p3,i)t, (^3,2)*, (P4,i)t 5 and {Pb)t- 
(P3,i)t = (cos(2po,i)( 3 Pi uw; + p\w Z - p\w 2 + 3p\uw\ - ip\uiw - 3p\w 2 w 1 

— 24piu 2 w + lObpiuw 3 + 6piuw 2 — §p\U\W\ + 6piu 2 w + 39p\W 5 

— 27piw 2 w 2 - 96pi4Pi i2 ww 2 - 24pi ! ipi, 2 u/ t + ^&P\,\P\.2WW 2 

— 2Api^pi^w\ + 36p3uw + 12p3U> 3 — 12p 3 w 2 — 12u 2 wi 
+ A8uu\w + 39uw 2 w\ + 3uiw 3 — 6uiw 2 + 6u 2 wi 

— 6u 3 w ~ 9w A W\ — 12w 2 Wz — \SwwiW 2 + 6u> 3 ) 

+ 2sin(2p ,i)(36pipi ; ipi i2 uw + 12p 1 p li ip li2 w 3 - 12p 1 p 1A p lt2 w 2 
+ 5Apiuwwi + 54piwiui 2 + 54piw 3 wi — I8P1WW3 + 12pi^pi^ 2 uwi 

— 12pi ; ipi i 2UiW — I2pi^ipi^ 2 w 2 wi — 9u 2 w 2 + 18uw 4 + 27uww 2 — 9uw\ 
+ 9u\ww\ + 3u 2 w 2 + 2w 6 — llw 3 ii;2 + 12w 2 wl — 3ww4 + 3wiw 3 — 3w 2 ) 
+ 6pi,2(— 3p\uw — p\w 3 + p\w 2 — 2p\uw\ + 2piuiw + 2piw 2 w\ 

+ 8u 2 w + uw 3 — 2uw 2 + 2uiW\ — 2u 2 w — w 5 — 3w 2 w 2 ))/12, 
(P3,2)t = (2cos(2po,i)(36pipi ! ipi !2 uw + l2pipi A pi^w 3 - V2p x p\,\P\fiWi 

— 54piuwwi — 54piuiw 2 — 54piw 3 w\ + 18piuw 3 + \2p\^p\^ 2 uw\ 

— \2p\^p\_ 2 u\w — 12pi^pi t2 w 2 w\ + 9u 2 w 2 — I8uw 4 — 27uww 2 + 9uw\ 

— 9u\wwi — 3u 2 w — 2w + llw w 2 — 12w w 1 + 3ww4 — 3wiw 3 + 3w 2 ) 
+ sin(2po,i)(3piMW + p\w — p\w 2 + Sp^wi — 3p 1 uiw — 3p 1 w w\ 

— 2Apxu 2 w + 105piuu> 3 + 6p\uw 2 — 6piUiWi + 6piu 2 w + 39piw 5 

— 27piw 2 w 2 + 96pi^pi, 2 uw 2 + 24pi ! ip^ 2 w 4 - 48pi^pi^ 2 ww 2 

+ 2Ap\ t ipi t2 w\ + 3&P3UW + 12pzw 3 — 12p 3 w 2 — 12u 2 wi + 48uuiw 
+ 39uw 2 wi + 3wiui 3 — 6u\w 2 + 6u 2 wi — 6u 3 w — 9w 4 wi — 12w 2 w 3 

— 18wwiw 2 + 6u> 3 ) + 6pi ! i(— 3p\uw — p\w 3 + p\w 2 — 2p\uw\ + 2pimw 
+ 2p\w 2 w\ + 8u 2 w + uw 3 — 2uw 2 + 2uiWi — 2u 2 w — w 5 — 3w 2 w 2 ))/12, 
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(P4,i)t = (8cos(2p ,i)(3p?Pi,2Mw+p^i,2W 3 -pIpi,2W 2 + 3piPi, 2 MWl 

— 3p(pi,2UiW - 3p\pi y2 w 2 wi + 36pipl A pi, 2 uw + 12pipl 1 pi t2 w 3 

— 12pipf ;1 pi ;2 W2 - 24pipi j2 w 2 w - 3pipi !2 ww 3 + 6piPi, 2 uw 2 - 6/>ipi i2 UiWi 
+ 6piPi,2U 2 w + Spipi^w 5 + $PiPi,2W 2 W2 + 12pl 1 pi t2 uwi - Ylp\ x px^u\w 

— YlpX yPxtfuPwx - 48p M p 2 2 uw 2 - 12pi,ip 2 ;2 w 4 + 2Ap lyl pl 2 ww 2 

— 12pi ; ip 2 2«^i + 72p\,\u 2 w 2 + 18p\ y \uw A — 54pi j iuww 2 + lSpi^uw 2 

— 18pi ; iUiW«;i — 24pi j iM 2 w 2 — 4pi j iw 6 — 32pi4U> 3 w 2 — 2Ap\,\W 2 w\ 
+ 6p\^ww4 — 6pi,iwiw 3 + Qpi y iwl + 3Qpi y2 p 3 uw + 12pi y2 p 3 w 3 

— 12pi j2 p 3 w 2 — 12pi j2 u 2 wi + A8pi y2 uuiw + 39p\ y2 uw 2 wi + 3p\ y2 uiw 3 

— 6pi y2 u 1 w 2 + 6pi, 2 u 2 wi - 6pi t2 u 3 w - 9pi i2 w 4 wi - 12pi y2 w 2 w 3 

— 18pi y2 wwiw 2 + 6pi y2 wl) 

+ 8sin(2p 0A )(3plpi tl uw + pfp lt iw 3 - p\pi y iw 2 + 3p\pi yl uw 1 - 3pjpi yl uiw 

— 3p\pi yl w 2 w 1 + 36pipi t ipl 2 uw + 12pipi ; ipf 2 w 3 - 12pipi ! ip 2 2 w 2 

— 24pipi ! iu 2 w — Spipi^uu! 3 + 6pipi4uw 2 — 6pipi i iUi«;i + Qp\P\ y \u 2 w 

+ 3piPi,iw 5 + §PiPi,iw 2 w 2 + ^8p\ yl pi y2 uw 2 + \2p\ 1 px y2 w 4 - 2Ap\ 1 pi y2 ww 2 
+ I2p 2 1:1 pi y2 w\ + 12p 1A pl 2 uw 1 ~ 12p 1A pl 2 u 1 w - 12p 1A p 2 12 w 2 W! 
+ 36pi y ipzuw + 12pi y ip 3 w 3 — 12pi y ip 3 w 2 — I2pi y iu 2 wi + 48pi y iUUiW 
+ 39p\ y \uw 2 wi + 3p\ t \Uiw 3 — 6pi ! iiiiw 2 + &p\,\u 2 wi — &pi t \u 3 w 

— 9p\ y \w A w\ — 12pi ! iw 2 u'3 — 18pi t iww\w 2 + &p\ t \w\ — 72pi y2 u 2 w 2 

— 18pi y2 uw 4 + 54pi y2 uww 2 — I8pi y2 uw 2 + 18p\ y2 uiwwi + 24pi y2 u 2 w 2 
+ ^Pi,2 w& + 32pi i2 w; 3 uj 2 + 2Api y2 w 2 w\ - 6pi y2 wwi 

+ 6pi,2Wiw 3 — 6pi t2 w 2 ) - 3p\uw — p\w 3 + p\w 2 - Ap\uw\ + <±p\u\w 
+ Aplw 2 w 1 - 72p\p\ l uw - 2Ap\p\ l w 3 + 2Ap 2 l p\ l w 2 - 12p\p\ 2 uw 

— 2Ap\p\ 2 w 3 + 2Ap\p\ 2 w 2 + 48p 2 u 2 w + Qp\uw 3 ~ \2p\uw 2 + 12plmwi 

— 12p\u 2 w — &p\w b — 18p\w 2 w 2 — 48pip 2 yl uwi + 48pip 2 yl uiw 

+ 48pip 2 1 w 2 w;i — 48pip 2 2 uwi + 48pip 2 2 uiw + 48pip 2 2 w 2 wi — IMpip^uw 

— 48piP3W 3 + 4:8pipzw 2 + 48piu 2 wi — 192piuuiw — I56piuw 2 wi 

— 12piuiw 3 + 24piuiw 2 — 24piu 2 wi + 24piu 3 w + 36piw 4 wi + 48piw 2 W3 
+ 72p 1 ww 1 w 2 - 24piw 3 + 192p 2 s u 2 w + 24p 2 11 uw 3 - A8pj A uw 2 

+ A8pl 1 u 1 w 1 - 48p 2 1 u 2 u; - 24p 2 1 «; 5 - 72p 2 1 w 2 w 2 + 192pl 2 u 2 w 
+ 2Ap\ a uw 3 - 48p 2 . 2 ww 2 + 48p 2 . 2 wiwi - 48p 2 . 2 u 2 w - 24p 2 2 w 5 

— 72p\ 2 w 2 w 2 — A8pzuw\ + 48p3Miw + 48p3Ui 2 wi — 252u 3 w — 36u 2 w 3 
+ 60u 2 w 2 - lMuuiivi + 240uu 2 w - 7uw 5 + 342uw 2 w 2 + 24uwwl 

+ 144u 2 w + 228uiw 2 wi + 94u 2 w 3 - 24u 2 w 2 + 24m 3 wi - 24u 4 u; 

— 13u> 7 + 47w 4 w 2 — 4w 3 wf — 78ui 2 u>4 — 6OWW1W3 
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- 120ww 2 + 60w 2 w 2 )/48, 

(p 5 ) t = (54m 4 + 180u 3 «7 2 - 12u 2 u 2 + 126u 2 w 4 - 282u 2 ww 2 - 12u 2 w 2 

— 84uuiivwi — 174uu 2 u> 2 + 6uu 4 + I8uw 6 — 300uw 3 w 2 — 90uw 2 w 2 
+ 66UWW4 + 61LW1W3 + 48uw 2 + A2u\w 2 — 6uiu 3 + \2u\w z w\ 

+ 42uiww 3 — 48uiWiW2 + 6u 2 — 39w 2 u> 4 + I62u 2 ww 2 — 48it 2 u> 2 
+ 48u 3 ww 1 + 2lu 4 w 2 - A2w 5 w 2 - 12w; 4 u7 2 + 35w 3 w 4 + 120w 2 wiw 3 
+ 195«; 2 w;2 + 120uw 2 w 2 — 6wwq — 30w 4 + 6wiw 5 — 6w 2 W4)/6. 
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